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Abstract 

The propagation of neutrinos in a gravitational field is studied. A method of 
calculating a covariant quantum-mechanical phase in a curved space-time is 
presented. The result is used to calculate gravitational effects on the neutrino 
oscillation in the presence of a gravitational field. We restrict our discussion 
to the case of the Schwarzschild metric. Specifically, the cases of the radial 
propagation and the non-radial propagation are considered. A possible ap- 
plication to gravitational lensing of neutrinos is also suggested. 
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I. INTRODUCTION 



Neutrino oscillations in a flat space-time have been extensively studied in the past by 
using both plane waves |L| and wave packets @ to represent the emitted neutrinos. In 
particular, it has been shown || that the standard treatment of neutrino oscillations in 
the plane wave approximation is valid only for extremely relativistic neutrinos, whereas for 
a general case, the wave packet treatment is essential. In this paper, we discuss how the 
results in a flat space-time are modified in a curved space-time. That is, we calculate the 
quantum mechanical phase of neutrinos that are produced and propagate in a gravitational 
field. Our derivation of neutrino oscillation formula in a gravitational field will be based on 
the covariant form of the quantum phase that arises due to the assumed mixing of massive 
neutrinos P|. First we consider the case of neutrinos that are emitted and propagate in 
a radial trajectory in the Schwarzschild metric. Such a gravitational effect can, in princi- 
ple, modify the standard vacuum oscillation formula for the solar and supernova neutrinos. 
Although the size of the effect is far beyond the current experimental detectability, in par- 
ticular for the solar neutrinos, it may certainly be of interest for the neutrinos from very 
massive sources. ( It is well known that gravitational influence on the MSW effect for the 
solar neutrinos is significant if the equivalence principle is violated ||.) In our derivation 
we have not assumed the weak field approximation. We then compare our results with the 
previous results in the case of the radial trajectory obtained by [§-0] with clarifying remarks 
on the differences in the results and interpretations. As a further application, we also con- 
sider the case of non-radial propagation in a gravitational field. Finally, we discuss possible 
gravitational lensing effects on the neutrino oscillations, for which it is necessary to resort 
to the weak field approximation. It is to be noted that in the last two cases, due to the 
angular spread of neutrinos with different masses in the presence of the gravitational field, 
a proper way to treat the neutrinos is to resort to the wave packet formalism. (However, 
a complete, covariant description of neutrino propagation in a gravitational field in terms 
of wave packets is beyond the scope of this paper. This issue will be addressed elsewhere 
||.) In order to discuss the problem in a transparent way, and in order to compare with 
the previous analyses, we restrict ourselves to the discussion of relativistic neutrinos, where 
a plane wave analysis can be employed. 

The plan of the paper is as follows. In Section II, we briefly review the standard treatment 
of neutrino oscillations in a flat space-time using the plane wave formalism. In Section III 
we extend the plane wave analysis to the case in the presence of a gravitational field. For 
definiteness, we discuss the neutrino oscillations in a field described by the Schwarzschild 
metric. Specifically the cases of radial and non-radial propagation of neutrinos are discussed. 
In the last part of Section III, we suggest the possibility of gravitational lensing of neutrinos 
and evaluate the resulting flavor-changing oscillation probability. 

II. NEUTRINO PROPAGATION IN A FLAT SPACE TIME 

Let us consider a neutrino produced at a space-time point A(tA,XA)- Since it is produced 
by a weak interaction process, it emerges as a flavor eigenstate \v a ), which is a superposition 
of the mass eigenstates \vk), i-e. 
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k) = E^k>> (!) 

where £7 is the unitary mixing matrix of the neutrino fields. What actually propagates are 
the mass eigenstates, whose energy and momentum are Ek and pk, respectively, and they 
are related by the mass-shell condition as 

El=p k 2 + m 2 k . (2) 

Both Ek and pk are determined by the energy and momentum conservation at the production 
point A and, in general, they are different for different mass eigenstates. In a flat space-time, 
the propagation of the state \uk) is described by a plane wave 

\v k (t,x)) = exp(-i$ k )\v k ) , (3) 

where 

$ fc = E k t - p k ■ x . (4) 

Neutrino oscillations take place due to the fact that different states \v k ) propagate dif- 
ferently because they have different energies, momenta and masses. When they arrive at 
a detector located at a space-time point B(t B ,x B ) which detects flavor eigenstates via a 
weak interaction process, they have developed a relative shift in their phases. In order 
for the oscillation to occur and to be observed, some requirements must be met. First, in 
addition to the standard assumption of mixing of massive neutrinos, the mass eigenstates 
must be produced coherently. This implies that the interference is possible only among mass 
eigenstates produced in the same process, because neutrinos produced by different processes 
have, in general, random relative phases in their wave functions, which destroy the coher- 
ence. Secondly, the states have to be detected at the same time t B and at the same place 
x B - 

Under these circumstances, the interference can take place and the oscillation phe- 
nomenon arises. The probability that the neutrino produced as \u e ) is detected as \u^) 
is, therefore, (in the case of two generations, where U is parameterized as a function of the 
mixing angle 9 in the usual way) || 

V[y e -> z/ M ) = |(z/> e (t B ,f B ))| 2 = sin 2 (2#)sin 2 (^j , (5) 

where $12 = $1 — $2 and $k (k = 1, 2) are the phases 



rtB r%B 

E k (t B - t A ) - p k ■ (x B - x A ) = E k dt-pk-l dx , (6) 



acquired by the mass eigenstates. 

The expression for the phase in Eq.@ can be written in a covariant form, which is 



suitable for the subsequent application in a curved space-time, as |10] 



$ k = [%l k) dx^ (7) 
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where 

pf = m ^^r , (8) 

is the canonical conjugate momentum to the coordinates x M and g^ u and ds are the metric 
tensor and the line element, respectively. This covariant phase in Eq.(|7]) was first discussed 
by Stodolsky [jlOfl , and has been used in |5|-[7j to calculate the neutrino oscillation phase 
difference. 

Equation @ represents the oscillation probability for a neutrino produced at the space- 
time point A(t A , x A ) and detected at a given space-time position B(tB,XB)- In actual 
experiments, however, the time difference (t# — t A ) is not measured, whereas the relative 
position \xb — %a\ of the source and the detector is known. In the plane wave formalism, 
this can be taken care of consistently only for relativistic neutrinos by replacing (t B — t A ) 
with UTT 



{tB ~t A ) ^ \x B ~x A \ , (9) 

thus the time difference does not appear in the formula for the oscillation probability. In 
this approximation, the phase shift of Eq. becomes 

<S>k = {E k - \p k \)\x B - x A \ . (10) 

Applying the relativistic expansion m k <C E k , we can approximate, to the first order, 



where Eq is the energy for a massless neutrino. Therefore, we have 

mi 



E k -\p k \ = E k -^El-ml^^-, (12) 



which leads to the standard result for the phase 
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®k-^\x B -x A \. (13) 
it,® 

The phase difference responsible for the oscillation can be given by Eq.(|i~3"D as 

where Am^ = m\ — m|. 

For more general situations, where some or all of the states v k are non-relativistic, the 
above discussion cannot be applied, and a wave packet analysis is required @. In this 
case, the relation in Eq.([|) is no longer valid, and moreover the problem of the coherence 
of the different states at the detection position has to be taken into account. However, for 
relativistic neutrinos, the wave packet formalism shows that the approximation of Eq.(|9]) 
is indeed appropriate, and the oscillation probability V{y e — > v^) has the form of Eq.(|), 
where the phase shift Q k j is given by Eq. (|I4"D . 
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III. NEUTRINO PROPAGATION IN THE SCHWARZSCHILD METRIC 



Let us now turn to the discussion of the propagation of neutrinos in a gravitational 
field. For the sake of definiteness, and also because it may represent a situation of possible 
physical interest, we will discuss the propagation in a gravitational field of a non-rotating 
spherically symmetric object, which is described by the Schwarzschild metric. The situation 
under consideration can be described by the line element in the coordinate frame {t, r, </>} 
as 

ds 2 = B(r)dt 2 - Bir^dr 2 - r 2 d$ 2 - r 2 sin 2 M$ , (15) 

where 

, s / 2GM\ . . 

B(r) = (l- _) , (16) 

and G is the Newtonian constant and M denotes the mass of the source of the gravitational 
field. Since the gravitational field is isotropic, the classical orbit may be confined to a plane. 
Hence, we can choose it to be on the equatorial plane = tt/2, and we have dd = 0. 
The relevant components of the canonical momentum of Eq. ® are 

(17) 
(19) 

and they are related to each other and to the mass m k by the mass-shell relation 

m\ = sTpMpP (20) 

1 ( C0\2 

= ^{pf ] f-B{r){ P ^) 2 -^. 
B[r) r z 

The fact that the metric tensor components do not depend on the coordinates t and 
(f) ensures that their canonical momenta p\ and p£ are constant along the trajectory. 

We define the constant of motion to be = pj> and Jk = —p^- They represent the 
energy and the angular momentum which an observer, located at r = oo, sees for the mass 
eigenstate v^. They differ from the energy and the angular momentum measured by an 
observer at a position or those at production point r^. The correct way to define the 
energies which are actually involved in a realistic situation is not, in general, unique. For 
example, for a neutrino produced in the almost stationary shock wave of a supernova, a local 
static reference frame for the production point seems appropriate. On the contrary, for 
neutrinos produced in the accretion disk around a black hole, a free-falling orbiting system 
seems proper. Similar arguments apply to the detection of neutrinos. For example, in the 
case of solar neutrinos, the detectors are in the free-falling frame. The general situation can 
be rather complicated and every case must be carefully dealt with. In our discussion we will 





= m k B(r) 


dt 
ds 




= -m k B 


(r)" 


pf 




dcf) 


= — rrikr 2 


ds 
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choose the local reference frame. The local energy, defined as the energy measured by an 
observer at rest at a position r, is related to E k from the transformation law which relates 
the local reference frame {x a } = {t, r , ip, 6} to the frame = {t, r, ip, 8} [|T!J 



x a — L a ll x^ 1 ; — L a ^ u ri & p , (21) 



where L a M are the coefficients of the transformation between the two bases: 



L t t = y\gtt\, L r r = ^J\g rr \, L# = ^\g m \, = ylflWl > others = . (22) 

Therefore, the local energy is 

Et C \r) = \9u\- 1/2 E k = B(r)-^E k . (23) 

In order to obtain the neutrino oscillation probability in a gravitational field, we will 
calculate the interference of the wave functions of different mass eigenstates created at a 
space—time point A and detected at a space-time point B. In the plane wave approximation, 
the phase of each mass eigenstate v k is defined by the covariant expression in Eq.(|7]) and 
the interference of the k th and j th mass eigenstates is given by the phase difference 

$ kj = f A {vf - vf) dx* = $ k - . (24) 

Here the integration must be made on a definite space-time trajectory from A to B. Follow- 
ing the standard treatment of the oscillations of the relativistic neutrinos in a flat space-time, 
as discussed in Section ||, we will calculate the interference phase in Eq. ([24]) along the light- 
ray trajectory from A to B. This corresponds to the approximation in Eq.(|9|) for the flat 
space-time case. We emphasize that the phases in Eq.(^4j) are not the phases on the classical 
trajectory of the mass eigenstates [|l(| but the phases calculated on the light-ray trajectory. 



We will see that for relativistic neutrinos the result for the phase difference in Eq. (|24|) is 
proportional to Am1j/2E, as in the standard treatment of neutrino oscillations in a flat 
space-time. 

We will now define the phase acquired by the mass eigenstate v k when it travels from 
point A(tA, ta-, 4>a) to point B(tB,rB, 4>b) as 

= [ B [E k dt - p k (r)dr - J k d(j>] , (25) 

J A 

where we have defined p k (r) = —p^ k \ The integration in Eq.(P5|) is performed along the 
light-ray trajectory which links the space-time points A and B. At this stage, we note that 
E k and J k , which are constants of motion for the geodesic trajectory of the k th eigenstate, 
are no longer constant along the light-ray trajectory. Instead, the energy at infinity E and 
the angular momentum J at infinity for a massless particle are constant along the light-ray 
path. Therefore, E k and J k cannot be taken out of the integration in Eq.(^) and some 
caution is required for the calculation. We will show explicitly in the following Subsections, 
however, that in the relativistic limit, this problem can be circumvented. 

We now discuss two different situations: radial propagation and non-radial propagation. 
In the last Subsection we will address the possibility of gravitational lensing. 
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A. Radial propagation 



For neutrinos propagating in a radial direction, we have d0 = and no angular momen- 
tum. Equation (|25|) is reduced to 



I' A 







dr , (26) 



dr 

where Pk(r) is obtained, from the mass-shell relation Eq.(^) with J k = 0, as 

1 



p k (r) = ±— - ^El - B(r)ml , (27) 



and the light-ray differential (dt/dr)o is 

( dt\ 



K dr J Q £>(r) 



(28) 



In Eqs.(^3) and (|28|) , the sign (±) apply to neutrinos propagating outward (+) or inward 
(— ) of the gravitational well, respectively. Therefore, the quantum mechanical phase is 

^ = ± / J - - 5(r)mg) ^ . (29) 

At this point, we apply the relativistic expansion using the energy at infinity E k as a reference 
value, i.e. m k -C -Efc. As in the flat space-time case, the following relation holds 

Ek ~ E ° + ° (if) ' (30) 

where Eq is the energy at infinity for a massless particle. Taking into account that < 
B(r) < 1, we have 



El - B(r)ml * E k - B{r) £l . (31) 

Then, the phase in Eq. (|29"D is approximated by 

$ fc ^±/ ^dr. (32) 
Jr A * J-<a 

Since the integration is performed along the light-ray trajectory, E Q is constant and the 
integration is easily performed to give 



2E 

The phase shift which determines the oscillation is, therefore, 



r B -r A \ . (33) 
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kj 



H 

2E 



r B - r A \ 



(34) 



We note that the derivation of this result does not depend on the weak-field approximation. 

The result for the phase shift in Eq.(K3J) is in agreement with that in M. The authors 
in [0] calculated the phase of each particle along its classical trajectory and then introduced 
an initial time difference in the phases for the states and Uj. Although they arrived at 
the correct result (Eq.(Q)), their derivation is not justified because the neutrinos produced 
at different times do not interfere because they possess initial relative random phases. We 
believe that the correct approach is to consider the interference between mass eigenstates 
produced at the same space-time position and detected at the same space-time point, related 
by the light-ray relation of Eq. fl28|) . On the other hand, any comparison of our result with 
that in || is problematic since the energy E used in || is not clearly defined among others. 

Some comments on the definition of "relativistic" neutrinos are in order here. Let us 
consider the following cases : 



(1) ml^El 



(2) m2« 

(3) m2« 



E 



E 



(loc) 

k 1 

(loc) 



r B 



relativistic at infinity, 
relativistic at the source, 
relativistic at the detector. 



In case (1), the ratio of m\ to any local energy E^ oc \r) is, from Eq. (|23"D and B(r) < 1, 



mi 



m 2 k 



m 2 k 



Et\r) 



-JiB(r) < — < 1 



(35) 



so that the neutrinos are even more relativistic at r < oo, and the approximation in Eq.(pl|) 
is certainly justified. 

Case (2) needs a caution when the observer happens to be at infinity, because the ratio 
of ml to the energy at r = oo becomes 



El 



mi 



B{r A ) 



(36) 



That is, even if neutrinos are produced highly relativistically, they are not guaranteed to be 



relativistic at r = oo, unless r A 3> 2GM 



1 



{m k /E k 



(loc) 



r A )f 



For the obvious reason 



that non-relativistic neutrinos cannot be detected, at least with known techniques (assuming 
that neutrino masses are much smaller than 1 MeV), however, the lack of validity of the 
relativistic condition at infinity not only means that the approximate formula ( |3i~D is not 
valid, but also that in practice such neutrinos are not detectable at infinity. 

Case (3) deals with an observer under the influence of a sizeable gravitational field. In 
this case, neutrinos stay always relativistic along their path {r A < r < r B ), which validates 
approximation in Eq.([U]), for we have 



B(r) 



mt 



B(r B ) 



Et\r B ) 



< 



mi 



Et C \r B ) 



< 1 



(37) 
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In short, neutrinos are assumed to be "relativistic" when they are relativistic at in- 
finity, relativistic at the detector, or relativistic at the production point with r A ^> 



2GM 



(m k /Et\r A )) 



and then Eq.(|5BD provides the correct quantum phase. 
As a final comment, we wish to compare Eq.(^) with that of the flat space-time case. 

) appear 



As they stand, the expressions of the phase in Eq.(33) and the phase shift in Eq.([ 
identical to those of the flat space-time case. However, the gravitational effects are present 
implicitly in Eq.(|33"D and Eq. (|34|) . In the absence of a gravitational field, E is the energy of 
the neutrino as seen by any observer along its trajectory, and (r B — t a ) is the distance over 
which a neutrino propagates. Therefore, Eq. (|33|) gives the standard result shown in Eq. (|13|) . 
However, in the presence of gravity, the propagation of a neutrino is over its proper distance 



L p = 



TD 



' dr 



(38) 




In Jr A - 2GM + Jr~ A 



+ 2GM 



To simplify the following discussion, we consider the case of a weak field, where L p is 
approximated to 



L ~ 



Tb - r A + GM In 



(39) 



This shows that, in a gravitational field, the effective length in the phase (i.e. (tb — r A )) is 
shorter than L p . Moreover, the energy measured by a detector at tb is not Eq, but rather the 
local value E^ loc \rB)- When expressed in terms of the local energy and the proper distance, 
the phase shift of Eq.([34|) in the weak field approximation is 



' Am lj L P 
2E^ oc )(r B ) 



1 - GM 



ln^- 

r A 



(40) 



The first parenthesis on the right-hand side in Eq.([40D is analogous to the flat space-time 
oscillation phase. The second parenthesis represents the correction due to the gravitational 
effects. 

The proper oscillation length L^f, in the weak field approximation, is 



T OSC 



{tb) 



ATiE Xoc (r B ) 



GM 



In 1 



4irE loc (r B )\ 47rE loc (r B ) 



Am 2 Tb 



+ 



Am 2 Tb 



(41) 



where the quantity in square brackets is negative. We conclude, therefore, that the proper 
oscillation length is increased in the gravitational field, as expected. 



B. Non radial propagation 

In this Subsection, we discuss the case of the propagation along a general trajectory. In 
contrast to the radial case, the motion has an additional angular dependence. The phase 
$fc is 



9 



TA 



E, 



'df 
dr 



~Pk{r) - Jk 



dr 



dr 



(42) 



where the integral is taken along the light-ray trajectory that links the production point 
A to the detection point B. In Eq . ((42|) , the quantities (dt/dr)o and (d<p/dr)o along the 
light-ray trajectory are 



E 



'dt\ 

K drJ Q B 2 (r)p (r) ' 

d(p\ J 1 

dr J n r 2 B(r)p (r) 



(43) 



It is convenient to express the angular momentum Jk as a function of the energy E k , the 

(44) 



impact parameter b and the velocity at infinity t>jf°' > |13| 



Since at r = oo the metric is Minkowskian (no gravity), we can write 

1 



^ 2 - ml 



Ek 



2EI ' 



(45) 



where in the last equality we used the relativistic approximation up to the order 0(m|/ 'El). 
The angular momentum of a massless particle, J , is obviously 

J = E b. (46) 

With Eqs.(|4l3|)-(f46"|), the expression of $^ in Eq.(|42]) can be conveniently arranged as follows 



dr- 



En 



lr A B(r)p (r) 
The mass-shell condition Eq.(f 



E k B(r)p (r) 



B(r) 
) gives 



En 



-Pk{r) 



E k b 2 



m k 



(47) 



B(r)p (r) = ±E Jl-B(r)^ 



B{r)p k {r) = ±E k 



\ 



b 2 



mi 



b 2 ' 



1 - B(r)— - B(r)-S 1 - 



El 



(48) 
(49) 



±E k \ l-B(r) 



b 2 



1 - 



B(r)(l-b 2 /r 2 ) ml 
l-B(r)b 2 /r 2 2E\ 



The last approximate equality in Eq.(f49|) is due to the relativistic expansion. In Eqs.(|48"D and 
(fl9l), the sign (±) is determined by whether dr is positive (+) or negative (— ). Substitution 
of Eqs.([|^) and Q49| ) into Eq.(f4"?D simplifies the expression for the phase Q k to 



/ 

J r 



~ ' dr 



En 



Ei 



ml 



r A B(r)p (r) K 2E 2 



(50) 
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Since the following relation holds, in the relativistic approximation of Eq. fl30|) , 



<3>fc can be expressed as 



rat 



m k 



Er 



2E t 



m k 
2E 2 ' 



dr 



J rA 



1 - B(r)(b 2 /r 2 ) 



(51) 



(52) 



Equation (|52j) is the phase acquired by the mass eigenstate \v^) for a non-radial propa- 
gation from the source A to the detector B. In the limit b ^ 0, which reduces the motion 
to be radial, Eq.(|33|) is recovered. We also notice that the integrand in Eq.(|52D is divergent 
at the point of the closest approach ro, defined by the condition that the rate of change of 
the coordinate r with respect to the angle <fi vanishes: 



dr 

- = 



El 



1 - - 2 B(r ) = 



(53) 



However, the integral which gives the phase is finite. We will show this explicitly in the 
weak field approximation. 

The expression of $^ obtained in Eq.(|52]) is valid for any spherically symmetric (and 
time-independent) field. It has been derived without any assumption on the strength of the 
gravitational field. In order to gain more physical insight, however, we perform a weak field 
approximation, which allows us to perform the integration analytically. The approximation 
is valid if the field is weak enough to satisfy the condition GM <C r for all the r's along the 
trajectory under consideration. For example, the gravitational field of the sun at its surface 
is about GM Q / R & ~ 2 x 10" 6 and that of a galaxy is about G(lO n M )/3Okpc ~ 1.6 x 1(T 7 , 
both of which justify the weak field approximation. Whenever the weak field approximation 
is applied, we keep the expansion up to the order 0(GM/r). 

First, let us consider the case where a neutrino is produced in a gravitational field and 
then propagates outward from the potential well non-radially. The weak field approximation 
allows us to expand 



1 - B(r 



b 2 



1- 



b 2 



1 + 



GM b 2 



b 2 



(54) 



The phase is then easily integrated and becomes 



mt 



2E n 



r B~b 2 -Jr 2 A 



b 2 + GM 




TA 



r\ — b 2 



(55) 



We notice again, as a consistency check, that the radial limit b — > gives the same expression 
as given in Eq . (|33|) . 

The second situation is when a neutrino moves around the massive object, crossing the 
closest approach point at r = r . Taking into account the sign of the momentum, the phase 
is 
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®k(r A -> r -> r B ) 



mi 



2E, 



r.4 



dr 



+ 



mi 



J ro 



yjl - B{r){b 2 /r 2 ) 2E o Jr ^1 - B{r){b 2 /r 2 ) 



TB 



dr 



(56) 



The position of the closest approach can be solved from Eq. (|53| ) in the weak field approxi- 
mation, as 



(57) 



Substituting Eq.(|57|) into Eq.(|5£j), we have 

,2 



2E 

ml 



r'i -r 2 n + \jr\ - r\ + GM 



lr A - r \r B - r 



rA + r 



rB + r 



(58) 



2E, 



o 



r\ — b 2 + \jr 



l-b 2 




+GM 



We observe that, in this case, the radial limit b — > is meaningless, because it would 
correspond to a radial motion which crosses the gravitational source, where our description 
becomes inadequate. 

For b <C ta,,b, Eq.(|58|) is reduced to (up to the order of {b 2 /r 2 AB )) 



2En 



J- a + r B ) 



1 - 



+ 



2GM 



2r A r B r A + r B 



(59) 



which will be used to discuss the gravitational lensing in the following Subsection. It is 
interesting to note that Eq.(|59"D has a gravitational effect which does not depend on the 
distance between the source and the detector (assuming that this distance is much larger 
than the impact parameter b). That is, this gravitational effect integrated along a trajectory 
which passes close to a gravitational center induces a constant phase shift 2GM(ml/2E Q ). 
Furthermore, this constant phase shift does not depend on how close the trajectory passes 
to the gravitational center. Therefore, if, for example, a neutrino travels relatively close 
to several well-separated gravitational centers, the net phase shift becomes the sum of the 
phase shifts induced by each gravitational center. 



C. Gravitational Lensing of Neutrinos 

Let us consider a gravitational lens which is located between a source and an observer 
but off the line connecting the two. A neutrino emitted from the source can travel along two 
different paths, the proper distances of which are different and give the quantum interference 
at the detector. 

Oscillations arise due to the interference not only between the mass eigenstates and Vj 
travelling along each path, but also between the mass eigenstates propagating along different 
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paths (for definiteriess, we denote them as long-path (L) and short-path (S) ). A neutrino 
produced as a flavor eigenstate |z/ e ) = cos#|i/i) + sin#|z/ 2 ) at the source A(t A , r A , <Pa), evolves 
into (we consider only two generations) 



\v e ,B) = N J2 {cos9ex P (-i^ ath )\v 1 ) + sm9exp(-i&% ath )\v 2 ) 

path=L,S 



(60) 



where N is the normalization constant. The flavor-changing oscillation probability at the 
detector is then given by 



V(v e -> v^) = \ (vp\v e , B)\ 2 
= - cos 2 9 sin 2 9 



(61) 



1 + cos($f - $f ) + 1 + cos($ 
{cos($£ - $f ) + cos($f - $f)} 
{cos($£ - $f) + cos($f - 



d>f) 



The phases &P ath in Eq.([H]) can be evaluated along the light-ray trajectories as shown in 
the previous Subsection. Substituting Eq . (|59|) into Eq.(^l|), we have 



V(v e 



sin 2 (29) 



Am 2 X 



2GM 



S6 2 



sin 



x cos 



4£ 

'm\X Ab 2 



cos 



+ sin 



4E 2r A r B ) 
, , Em 2 X Ab 2 \ 
V 4E 4r A r B J 



sm 



X 4r A r B J 
' m\X Ab 2 \ 
v AE 2r A r B ) 
' Am 2 X Ab 2 \ 
, 4E 4r A r B 



(62) 



b\-b\ 



where we have defined X = r A + r B , Am = m 2 — mf, Hm = m 2 + mf, Ab = u L — u s 
and S6 2 = b\ + b 2 s . In the symmetric case where the lens is aligned with the source and the 
detector, Ab 2 = and S6 2 = 2b 2 and the above flavor-changing probability is reduced to 
that of the non-radially propagating neutrinos 



(Af> 2 =0) 



sin 2 (29) sin 2 



Am 2 X 

4£n 



1 + 



2GM 



X 2r A r B J 



(63) 



which can be obtained directly from Eq.(^). This is expected, since the symmetric case is 
equivalent to the case of the non-radial propagation. 

Obviously, the proper way to discuss the gravitational lensing effects on the neutrino 
oscillations would require a wave packet formalism. Such a study is beyond the scope of the 
present paper and will be given elsewhere. 



IV. CONCLUSION 

We have studied the propagation of neutrinos in a curved space-time and the modification 
to the neutrino oscillation by calculating a covariant quantum mechanical phase The 
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gravitational field considered in this work is that of a non-rotating spherically symmetric 
object, described by the Schwarzschild metric. Furthermore, we have assumed that neutrinos 
are relativistic so that a plane wave analysis can be applied. 

Radial and non-radial propagation have been discussed in the light-ray approximation. 
Although our phase for the radial motion is in agreement with the result of the previous 
work 0, the interpretations are different. Any comparison of our result with that in || is 
problematic since the energy used in is not clearly defined. 

The calculated phase appears identical to that of the flat space-time case. This is because 
the phase is expressed in terms of the asymptotic energy E and the coordinate distance. 
However, the gravitational effects do appear in the leading order if we express the phase 
with the locally measured energy and proper distance. As in the radial case, the phase 
of relativistic neutrinos for the non-radial motion has been obtained without resort to the 
weak field approximation. Assuming that the gravitational field is weak enough and the 
source and the detector are at a sufficiently large distance from the massive object, the 
phase is reduced to a simpler form as given in Eq. fl59|) . Finally, we have considered the 
gravitational lensing of neutrinos, i.e. the quantum interference when neutrinos propagate 
through different paths, and have derived the flavor-changing probability V{y e Vy) as 
given in Eq.(^|). 

Even though the measurement of the gravitational effects on the propagation and os- 
cillations of neutrinos is not feasible at present, we think that the understanding of these 
behaviors themselves is of interest. 



Acknowledgements. NF gratefully acknowledges a fellowship from the Istituto 
Nazionale di Fisica Nucleare, Italy. 



Note added. After the completion of this paper, we became aware of the paper by C.Y. 
Cardall and G.M. Fuller fl4]], which discussed a similar subject and obtained the results for 



the radial motion similar to ours. 
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